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Precision cosmology and big-bang nucleosynthesis mildly lavor extra radiation in the universe be- 
yond photons and ordinary neutrinos, lending support to the existence of low-mass sterile neutrinos. 
We use the WMAP 7-year data, small-scale CMB observations from ACBAR, BICEP and QuAD, 
the SDSS 7th data release, and measurement of the Hubble parameter from HST observations to 
derive credible regions for the assumed common mass scale ttIs and effective number A^a of thermally 
excited sterile neutrino states. Our results are compatible with the existence of one or perhaps two 
sterile neutrinos, as suggested by LSND and MiniBooNE, if nia is in the sub-eV range. 

PACS numbers: 14.60.St, 14.60.Pq, 98.80.Es, 98.80.Cq 



Introduction. — Big-bang nucleosynthesis (BBN) has 
for many years provided the best inference of the ra- 
diation content of the universe. This inference imphes 
the existence of approximately three neutrino flavors and 
hence three generations of elementary particles P, Q ■ Of 
course, by the early 1990s, precision measurements of the 
Z° decay width at LEP had superseded this early cos- 
mological evidence. More recently, observations of the 
cosmic microwave background (CMB) anisotropics and 
the large-scale structure (LSS) distribution have allowed 
us to probe the radiation density at the CMB decoupling 
epoch. Since this epoch corresponds to a much later time 
than the epoch of nucleosynthesis, inference of the radi- 
ation content of the universe from precision CMB/LSS 
measurements provides an important overall consistency 
check. The results are surprising on several counts. 

Firstly, when analyzing cosmological data in extended 
scenarios that allow for extra radiation, one consistently 
finds a trend towards a radiation excess, and this hint 
has become stronger over the years ■ The cosmic ra- 
diation content is usually expressed in terms of iVoff , rep- 
resenting the effective number of thermally excited neu- 
trino species. The standard value iVofi = 3.046 slightly 
exceeds 3 because e"^e~ annihilation provides residual 
neutrino heating Most recently, the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) collaboration found 
a 68%-credible interval of N^s = 4.34 



-0.88 



based on 

their 7-year data release and additional LSS data Q. 
A similar study including the Sloan Digital Sky Survey 
(SDSS) data release 7 (DR7) halo power spectrum found 



eft 



4.78^1^1 at 95% confidence see also Ref. 



Most recently the Atacama Cosmology Telecope Collab- 
oration has arrived at similar conclusions 

Secondly, for a good part of the past two decades, BBN 
and the observed primordial "^Hc abundance suggested 
^off ^ 3 [l3- A large body of literature was devoted 
to particle-physics constraints, assuming that additional 



radiation required other novel ingredients such as a Vg- 
De asymmetry [l|, 01- However, the BBN situation has 
also changed of late: Two recent studies find a somewhat 
higher -^He abundance of Yp = 0.2565 ± 0.001 (stat) ± 
0.005 (syst) [HI and Yp = 0.2561±0.011 [13 respectively, 
suggesting additional radiation during the BBN epoch. 

Taking these hints seriously, one quickly finds that it 
is not easy to account for additional radiation because 
well-motivated candidates are lacking. The strong dilu- 
tion of any pre-existing radiation at the quark-hadron 
phase transition at T ~ 170 MeV suggests that a novel 
radiation component should be produced after the QCD 
epoch, although there are counter-examples [isj . If neu- 
trinos are Dirac particles there is no obvious way to equi- 
librate the right-handed states because one needs larger 
couplings (e.g., magnetic moments or right-handed cur- 
rents) than allowed by other evidence [lj|. Thermal ax- 
ions can provide only a fractional contribution to A^off 
[ist . In contrast to earlier findings, it was recently shown 
that standard neutrinos with primordial asymmetries can 
provide any value for N^s [l6[. No simple mechanism, 
however, produces the required asymmetries. 

If low-mass sterile neutrinos exist and mix with active 
flavors, they can be thermally excited by the interplay of 
oscillations and collisions [l7[. The prior theoretical mo- 
tivation for such states is perhaps not strong. However, 
the signal for flavor conversion in the LSND ex- 

periment 18| , if interpreted in terms of flavor oscillations, 
requires low- mass sterile states U21|. Later, the Mini- 



BooNE experiment observed excess events in the v^^ — s- v,, 
channel that were not compatible with LSND in a two- 
flavor oscillation model [H [ll] . Early MiniBooNE data 
in the — >■ Df, channel were statistically not significant 



enough to confirm or refute the LSND signature [24 1. 
However, based on 70% more data, a clear excess was 
recently reported psj . If interpreted in terms of flavor 
oscillations, these data require CP violation and thus a 
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minimal scenario would involve 2 sterile neutrinos that 
can interfere appropriately [26, 27|, although this inter- 
pretation causes significant tension with disappearance 
experiments, notably atmospheric neutrinos. An alter- 
native is 1 sterile state together with nonstandard inter- 
actions [2^ . Yet another motivation for low- mass sterile 
neutrinos arises from a combined fit of MiniBooNE and 
Gallium radioactive source experiments [2^ . 

These scenarios together with the cosmological hints 
for extra radiation motivate us to reconsider the impact 
of low- mass sterile neutrinos on cosmological observables. 
Some recent studies have looked at light sterile states in 
the context of precision cosmology 30l432l| , but with em- 



TABLE I; Priors and standard values for the cosmological fit 
parameters considered in this work. 



phasis on providing constraints on such models. Further- 
more, these papers were written at a time when precision 
cosmology provided only very weak evidence for N^g > 3. 

Scenarios. — Thermal excitation of sterile neutrinos is 
complicated because of possible resonant matter effects 
and CP violating phases. However, for the eV mass range 
and the relatively large mixing angles needed to explain 
the oscillation experiments, the sterile states are likely 
strongly thermally excited 31, ll^l- Note however that 



the presence of even a small neutrino asymmetry can 
lead to incomplete thermalization [s^, [s^ . 

In the following we take the number of thermally ex- 
cited sterile neutrinos, TVs, to be an adjustable cosmo- 
logical fit parameter, while the active neutrinos are as- 
sumed to have a fixed abundance {Ni, = 3.046) so that 
A^e£F = 3.046 -I- Ns- For the neutrino masses, we consider 
two schematic scenarios. In the first scenario, the ordi- 
nary neutrinos are taken to be massless, while the sterile 
states have a common mass scale rus which is free to vary 
(the "3-1- A^s" scenario). The second scenario, the 
scenario, consists of massless sterile states and active neu- 
trinos that have an adjustable common mass mj^, i.e., all 
active species are treated as degenerate in mass. This is 
a good approximation because current cosmological data 
are not sensitive to the small mass splittings. 

Cosmological analysis. — Besides the new parameters 
A"s and nis or , we use a cosmological parameter space 
consisting of the standard "vanilla" ACDM parameters: 
The baryon density Wb = ^^b^^, cold dark matter density 
Wcdm = ^cdmh'^j Hubble parameter Hq, scalar fluctuation 
amplitude Ag, scalar spectral index rig, and the optical 
depth to reionization t. We use either or mg as a 
fit parameter, depending on the scenario under consid- 
eration, from which wc calculate the contribution to the 
matter density as (i) ili,h^ ~ Ng x ms/(93 eV) in the 
3 + Ns scheme, and (ii) n^h'^ = 3.046 x m^/{93 cV) in 
the A's -|- 3 case. Wc impose flat priors on all parameters, 
as detailed in Table HI Parameter estimation is performed 
using a modified version of the CosmoMC package [36 1. 

We use CMB data from WMAP after 7 years of obser- 
vation 0, as weU as from the ACBAR [sg, BICEP 
and QuAD (soj experiments. In addition, we use the halo 
power spectrum extracted from the SDSS-DR7 luminous 
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Prior 



h 
T 



ln(10i°As) 

Us 

nis, nil, (eV) 



0.01-0.99 

0.005-0.1 

0.4-1.0 

0.01-0.8 

2.7-4.0 

0.5-1.5 

0-3 

0-7 



red galaxy sample |40|. Finally, we impose a prior on 
the Hubble parameter based on the Hubble Space Tele- 
scope observations (Zlj . Since recent type la supernova 
(SN) luminosity distance data are plagued by unresolved 
systematic issues associated with the light-curve fitting 
methods (4^ . we do not use them in our default analysis. 

In Figure [TJ we show 2D confidence contours for the 
3 4- A's and A's -|- 3 cases. We also provide ID confidence 
intervals for Ag and rriu or for different scenarios in 
Table mi For A's -|- 3 (m^ = and 7^ 0), the evidence 
for A"s > is beyond 95%, although it is weaker in the 
opposite scenario 3-1- A'g where the new states are massive. 

The evidence for A"s> increases if wc include the full 
SN sample from Ref. [H] using the MLCS2K2-inferred lu- 
minosity distances. However this effect is less pronounced 
if adopting instead the SALT-ii light-curve fitting method 
(Table HH) . Independently of SN data, excess radiation is 
preferred even by our default data sets. 

BBN. — The relatively large A"s values compatible with 
this analysis are somewhat restricted by the primordial 
light-element abundances. Wc summarize current BBN 
constraints on As in Table [1111 For the "^He abundance we 
consider both the value found by Izotov and Thuan ll|, 
Y^'^ = 0.2565 ± 0.001 (stat) ± 0.005 (syst) and the value 
derived by Aver et al. = 0.2561 ± 0.011. We 

adopt the deuterium abundance of Pettini et al. [ist . 
log(D/H)p = 4.56 ± 0.04 and the baryon density in- 
ferred from our CMB-f-SDSS-hHST analysis, uj^^^ = 
0.02239 ± 0.00048. Theoretical primordial element abun- 
dances are calculated using PArthENoPE [43 |. 

It is noteworthy that independently of ^Hc, deuterium 
alone also provides nontrivial limits when combined with 
^CMB ^ q£ these results likely contain significant sys- 
tematic errors. We therefore refrain from a common like- 
lihood analysis and simply remark that both ^He and D 
somewhat prefer As > 0, but depending on the assumed 



errors for Yp it may be difficult to accommodate As = 2. 
Of course the oscillation data do not necessarily imply 
two fully thermalized sterile neutrinos. Indeed, Ref. [3l| 
shows that the degree of sterile state thermalization in a 
3 active -I- 2 sterile scenario depends strongly on the neu- 
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TABLE II: ID marginalized bounds on A^a and neutrino 
masses. In rows 3-6 we have used As = 1 or 2 exactly. Two- 
tailed limits are minimal credible intervals. 



Scenario 



Range for As 
68% 95% 



Range for rris 



or rriu 
95% 
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0.39-2.21 


< 3.10 


0.01-0.34 eV 


< 0.66 cV 


i Vs o 
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0.05-3.75 


< 0.22 eV 


^ n 49 pV 


2+3 






< 0.20 eV 


< 0.30 eV 


3+2 






< 0.29 eV 


< 0.45 eV 


1+3 






< 0.16 eV 


< 0.24 eV 


3+1 






< 0.35 eV 


< 0.48 eV 


Including 


supernova 


data (mlcs2k2); 




3 + As 


1.24-3.36 


0.26-4.31 


0.17-0.47 eV 


0.09-0.64 eV 


Including 


supernova 


data (salt 


-II): 




3 + As 


0.02-1.54 


< 2.57 


< 0.28 eV 


< 0.66 eV 


TABLE III: BBN constraints on 


As, using As > as a prior. 


Maximum 


of the marginalised 


posterior and 


minimal 95% 


credible interval (C.I 
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Data 




Posterior max 


95% C.I. 


Fpi'r+(D/H)p 


0.68 




0.01-1.39 


Fp^+(D/H)p 


0.69 




< 2.42 


(D/H)p+cUb°^« 


0.49 




< 2.12 



FIG. 1: 2D marginalized 68%, 95% and 99% credible regions 
for the neutrino mass and thermally excited number of degrees 
of freedom As. Top: The 3 + As scheme, in which ordinary 
neutrinos have = 0, while sterile states have a common 
mass scale ttIs. Bottom: The As + 3 scheme, where the sterile 
states are taken to be massless ms = 0, and 3.046 species of 
ordinary neutrinos have a common mass rriv 



trino mixing and mass parameters. Another possibility 
is the presence of a small lepton as ym metry, which can 
reduce the thermalization efficiency 34 , 3^ . Yet another 
option is that the oscillation data are explained by 1 ster- 
ile state plus new interactions [2^ . Still, for reference we 
provide mass bounds in Table for the cases of A's = 1 
or 2 exactly, besides the variable N^- 

Discussion. — Allowing for extra radiation as a cos- 
mological fit parameter, current cosmological data favor 
additional radiation compatible with recent hints from 
BBN. Assuming ordinary neutrinos to have a common 
mass nil, and the extra radiation to be massless, evidence 
for A''s > exceeds 95%, whereas the most constraining 
upper limit comes from BBN. With currently favored *He 
and D abundances, it would be difficult to accommodate 
two fully thermalized additional neutrino states. 

The usual degeneracy between extra radiation and the 
ordinary neutrino mass (Figure [1} weakens the neutrino 
mass limits, with ID credible intervals given in Table HIl 

However, it is more interesting to assume essentially 
massless standard neutrinos and attribute a possible 



mass to sterile neutrinos (3 + scenario). If we assume 
A'^g = 1, the 95% allowed mass range is nis < 0.48 eV. For 
A^s = 2 it is 0.45 eV (Table HJ), ahhough this case would 
be disfavored by BBN. For A's < 1, the 2D marginalized 
posterior probability distribution has a long tail so that 
ms ^ 1 eV is marginally allowed: the fewer sterile states 
there are, the larger the mass they can possess. 

The relatively small masses favored by cosmology are 
not assured to provide good fits to the short-baseline ap- 
pearance experiments — in principle a combined analysis 
as in Rcf. !31] is desirable but complicated because of 
the many parameters involved. Moreover, the degree of 
thermalization of the additional states would have to be 
considered in detail. Two fully thermalized states, cor- 
responding to A's = 2, are difficult to accommodate in 
BBN even with the new helium abundances. 

Our main message is that on present evidence, cos- 
mology does not exclude sterile neutrinos if they are not 
too heavy and thus do not contribute excessive amounts 
of hot dark matter. Quite on the contrary, both BBN 
and precision observations would happily welcome some 
amount of additional radiation corresponding to around 
one new thermal degree of freedom. Low-mass sterile 
neutrinos are one natural possibility. 

Low-mass sterile neutrinos mixed with active ones can 
strongly modify the neutrino signal from a core-collapse 
SN and r-process nucleosynthesis in the neutrino-driven 
wind (m, lisl-lisi . These effects should be studied in the 
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presence of two sterile flavors and CP violating phases. 

The Planck spacecraft, currently taking CMB data, is 
expected to boost precision determination of cosmologi- 
cal parameters. It is foreseen that Planck data will con- 
strain the cosmic radiation content at CMB decoupling 
with a precision of AiVoff ~ ±0.26 or better, so extra ra- 
diation at the level of 2 extra species could be detected 
with high significance [i^, Conversely, if Planck 

should find A^off = 3 ± 0.26, this would provide strong 
arguments against sterile neutrinos. Thus, the ongoing 
Planck measurements will make or break the new-found 
friendship between cosmology and sterile neutrinos. 
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